A Study of Structure Formation and Reheating in the D3/D7 Brane Inflation Model 
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We study the spectrum of cosmological fluctuations in the D3/D7 brane inflationary universe with 
particular attention to the parametric excitation of entropy modes during the reheating stage. The 
same tachyonic instability which renders reheating in this model very rapid leads to an exponential 
growth of entropy fluctuations during the preheating stage which in turn may induce a large contri- 
bution to the large-scale curvature fluctuations. We take into account the effects of long wavelength 
quantum fluctuations in the matter fields. As part of this work, we perform an analytical analysis 
of the reheating process. We find that the initial stage of preheating proceeds by the tachyonic in- 
stability channel. An upper bound on the time it takes for the energy initially stored in the inflaton 
field to convert into fiuctuations is obtained by neglecting the local fluctuations produced during 
the period of tachyonic decay and analyzing the decay of the residual homogeneous fleld oscillations, 
which proceeds by parametric resonance. We show that in spite of the fact that the resonance is of 
narrow-band type, it is sufficiently efficient to rapidly convert most of the energy of the background 
fields into matter fiuctuations. 

PACS numbers: 98.80.Cq 



INTRODUCTION 

In recent years a lot of interest has focused on the in- 
terface area between superstring theory and cosmology. 
The reasons are twofold. Firstly, cosmology can provide a 
possible arena to test string theory observationally. Sec- 
ondly, superstring theory may be able to resolve the con- 
ceptual problems from which the current realizations of 
scalar field-driven inflation suffer (see e.g. [l| for a dis- 
cussion). In particular, a lot of attention has recently 
been devoted to attempts to obtain periods of inflation- 
ary expansion of space in compactifications of superstring 
theories to four space-time dimensions. Since such com- 
pactifications contain a large number of scalar fields and 
since some of those remain massless before supersymme- 
try breaking, string theory gives rise to promising candi- 
dates for the inflaton, the scalar field driving the infla- 
tionary expansion of space (see e.g. 0, |3, \A Q for recent 
reviews of attempts to obtain inflation from string theory 
compactifications) . 

D-branes have played a particularly important role in 
the recent approaches to obtaining inflation from string 
theory ^6|, lZ, ^8,, ^, 10. . 11.] . A particularly promising model 
is the D3/D7 brane inflation model, in which the inflaton 
is the separation between the two branes in the direc- 
tions transverse to both branes [l2| • This model features 
a shift symmetry [l^ for the inflaton field which ensures 
that at the classical level, the potential is flat along the 
inflaton direction (this shift symmetry is broken under 
quantum corrections [l^). Inflation is induced by super- 
symmetry breaking terms. Various aspects of this model 
have been studied previously 15, 16, 17]. 

The large number of light moduli fields in string infla- 
tion models fo^], however, leads to a new danger: light 
fields which are not the infiaton are potential entropy 



modes and can give rise to en trop y fluctuations. The 
"curvaton" [2ll 22 [21, [H, ^ and "modulated re- 
heating" 271, [23, |29|, |30|, [SlJ scenarios are examples of this 
effect. Non- parametric generation of fluctuations from a 
global symmetry breaking (also invol ving an key isocur- 
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vature field) has been investigated in 

As pointed out in (s^j HI] , these entropy modes may 
undergo parametric instability during the initial stages 
of reheating. In particular, this instability occurs on 
super-Hubble (but sub-horizon) scales. This is possible 
[36( since the background fields carry the causal informa- 
tion to scales larger than the Hubble radius. The en- 
tropy fluctuations, in turn, will seed a curvature mode 
(which we will call "secondary" mode). There are two 
field scalar field toy models in which this secondary cur- 
vature mode dominates over the primary mode, the pure 
adiabatic linear perturbation theory mode. In fact [btI ] 
the secondary mode can grow to be larger than unity 
before back-reaction shuts off the parametric instability 
which drives the growth of the entropy fiuctuation. In 
this case, the model is ruled out by the observed (small) 
amplitude of curvature fiuctuations on large cosmological 
scales. 

The D3/D7 brane infiation model is a prototypical ex- 
ample in which there is an entropy mode which can un- 
dergo resonance durin g th e initial stages of reheating, 
the preheating 38, 3^ iO, 41] phase. In this paper, we 
study the growth of this entropy mode and calculate the 
magnitude of the induced secondary curvature fluctua- 
tions. It turns out that for reasonable particle physics 
parameters, the secondary fluctuations remain in the lin- 
ear regime at the end of the reheating phase. However, 
their amplitude can be (in the absence of back-reaction 
effects) larger than the amplitude of the primary mode. 
This would necessitate changing the parameters of the 
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model in order to reproduce the observed magnitude of 
the large-scale cosmological perturbations. 

Note that a similar conclusion was recently found ji^l 
in another model of brane inflation, the KKLMMT model 
llj. In that model, it had already been observed that 
due to their large phase space enhancement, second or- 
der fluctuations may give a dominant contribution [i^. 
We should emphasize that the effect we are calculating in 
this paper is an effect which arises in linear cosmological 
perturbation theory. Both the linear effects described 
here and the quadratic processes discussed in [i^ and 
elsewhere (e.g. [3, liH) can be operational at the same 
time. In fact, the enhanced growth of the linear pertur- 
bations which is the focus of this work will increase the 
strength of the quadratic processes. We expect that ef- 
fects similar to those discussed in will also arise in 
the D3/D7 system. 

The structure of this paper is as follows: we first re- 
visit the calculation of the amplitude of the spectrum 
of cosmological perturbations in the D3 /D7 inflationary 
model based on the magnitude of the primary adiabatic 
mode alone. Next, we give an approximate but analytical 
study of reheating in the D3/D7 brane inflation scenario. 
The initial stages of reheating are governed by a tachy- 
onic decay channel. The same tachyonic resonance chan- 
nel leads to exponential growth of super-Hubble entropy 
modes, which then in turn induce a secondary curvature 
mode. The last major section of this paper focuses on 
the determination of the amplitude of this mode. 

Our analysis of the reheating phase provides results 
which are interesting in their own right. A new aspects 
of our analysis is that we allow for an initial offset in the 
value of the "waterfall" field. Such an offset could arise 
as a consequence of the back-reaction of long wavelength 
[99! fluctuations of the matter fields. Such an offset will 
suppress the formation of domains on small scales. If the 
offset were sufficiently large, it could prevent the onset 
of the tachyonic instability. However, if we take the off- 
set to be generated by the abovementioned fluctuations, 
we find that for most of the realistic space of parameters 
of the D3 /D7 inflationary model, the resulting values of 
the field lie in the configuration space region for which 
the tachyonic resonance channel [46|] is effective. Thus, 
the initial stages of the transfer of the energy in the in- 
flaton fleld to matter fluctuations will proceed via the 
tachyonic instability in a time interval which is less than 
the characteristic time it takes for the matter fields to 
oscillate about their minima. A substantial part of the 
initial inflaton energy is converted to localized nonlin- 
ear fluctuations during the period of tachyonic resonance. 
Neglecting the back-reaction of these fluctuations on the 
homogeneous background, we see that a fraction of order 
unity of the initial infiaton energy density still remains 
in the background when the tachyonic instability shuts 
off. We are interested in demonstrating analytically that 
the transfer of energy from the inflaton to fluctuations is 



effective in rapidly draining almost all of the energy from 
the homogeneous fields. The dominant energy conversion 
processes are those related to the interactions of the non- 
linear fluctuations. These have been studied numerically 
in great detail in |46| . Since our aim is to demonstrate an- 
alytically that within less than a Hubble time the energy 
can be effectively drained from the homogeneous back- 
ground, we will neglect the nonlinear interactions and 
focus on the evolution of the homogeneous component 
of the fields once the tachyonic resonance shuts off (due 
to our approximation scheme, we will thus only obtain 
an upper bound for the decay time). We show that the 
decay of the residual homogeneous component of the in- 
flaton field occurs via a parametric resonance instability 
[H, [3^ 0, 41 1 . In spite of the fact that the resonance 



turns out to be of narrow-band type, it leads to an ef- 
ficient conversion of most of the energy from the back- 
gromid inflaton fleld to matter fluctuations. 

One of the advantages of the D3/D7 brane infiation 
model is that the matter fields are directly associated 
with the branes whose separation provides the inflaton 
fleld [lO0l |. As a consequence, there is a direct coupfing 
between the infiaton field and the matter fields, rendering 
it easy, as we will see, to obtain efficient reheating. This 
is a feature the D3/D7 brane infiation model shares with 
the D4/D6 model in which reheating was analyzed in 
and with the brane-antibrane inflation models in which 
reheating was studied in (see also [13), but con- 

trasts with the situation in inflationary models in which 
the infiaton dynamics and the standard model fields live 
in different throats, in which case efficient reheating is 
more difficult to obtain (see e.g. [H, [H 0, l56|, l57| ) . 



The outline of this paper is as follows. In the following 
section, we review the D3/D7 brane inflationary model, 
with particular emphasis on the form of the potential 
for the inflaton fleld. In Section 3 we discuss the vac- 
uum manifold of the model, and in Section 4 we revisit 
the calculation of the spectrum of cosmological perturba- 
tions. Section 5 focuses on the dynamics of the inflaton 
and of the other relevant scalar fields of the model before 
the onset of reheating, introduces the conditions for effec- 
tiveness of tachyonic resonance and demonstrates that for 
most of the relevant parameter space of the theory, the 
tachyonic resonance condition is satisfied. However, the 
tachyonic instability turns off before most of the energy 
of the infiaton field is released as matter fluctuations. To 
obtain an upper bound on the time it takes for the infla- 
ton energy stored in the residual oscillations to further 
dissipate, we study the later dynamics of the homoge- 
neous field components which proceeds via the less effi- 
cient parametric resonance channel (Section 6). Section 
7 is devoted to the computation of the secondary curva- 
ture fluctuation mode, the mode induced by the entropy 
fluctuations. 
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THE D3/D7 BRANE INFLATION MODEL 

The simplest way to construct a D3/D7 inflation- 
ary model in type IIB theory is to place a D3 brane 
(which we take to be in the — 3 directions) parallel 
to a D7 brane (which is in the — 3,6 — 9 directions) 
in flat space-time and break the resulting Af = 2 su- 
persymmetry completely by switching on non-primitive 
gauge fluxes on the D7 brane. The BPS condition is bro- 
ken, and the D3 brane starts to move towards the D7 
brane. The distance between the D3 brane and the D7 
brane in the 4, 5 plane appears on the world- volumes of 
the branes as a complex scalar S, whose value changes 
as the distance between the branes changes. This scalar 
field is identified with the infiaton, and the motion of the 
D3 brane towards the D7 brane is the slow-roll of hybrid 
infiation. 

Once the D3 brane comes close to the D7 brane, the 
open string between the branes become tachyonic be- 
cause of the non-primitive fiuxes on the seven brane. 
This is the top of the hybrid inflationary scenario where 
one would expect a tachyonic phase transition. The D3 
brane then falls onto the D7 brane and dissolves as a 
non-commutative instanton [l2]. The non-commutative 
instanton is not point-like and therefore the final config- 
uration of a dissolved instanton is non-singular [58]. One 
can also argue that in the end J\f = 2 supersymmetry can 
be restored. 

This simple picture has many inherent problems. 
Firstly, being non-compact, the gravitational solution is 
ten dimensional and a direct dimensional reduction will 
give us vanishing Newton's constant in four dimensions. 
One way out is to compactify the other six dimensions. 
Secondly, once we compactify, new problems have to be 
faced. An immediate issue that on a compact space, 
charge neutrality must be enforced, which requires in- 
troducing additional elements to the construction. Next, 
the problem of moduli stabilization arises, because with- 
out stabilizing the moduli, the system will decompactify 
back to the ten dimensional solution via the so-called 
Dine-Seiberg runaway [59j. Thirdly, one needs to avoid 
an over-production of cosmic strings when the D3 dis- 
solves in the D7 brane. 

Ways to solve these problems were addressed in the pa- 
pers 13, 17 1 where the proposal was to consider a warped 
solution of the form 
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e^(^'*^5m«d2/"d2/" , (1) 



where A, B are the warp factors and are the coordi- 
nates of the internal space. In [l^ the metric gmn was 
chosen to be the metric of KZ x T"^ jZi with being 
a combination of an orbifold and an orientifold opera- 
tion (see 12| for details). Both the D3 and the D7 are 
points on the = T^/Z2, and the motion of the D3 is 
trigerred by three form fluxes H^^s and iJB>B |101j . The 



D3 brane is oriented along x^'^-'^''^ and therefore is also a 
point on the K3 manifold. The D7 brane wraps the K3 
and is oriented parallel to the D3 as we discussed before. 
Supersymmetry is broken either by making 



(2) 



where G = Hrr + tHmsi or by switching on ISD (1,2) or 
(0,3) forms (recall that ISD (1,2) fiuxes also break susy. 
For more detailed analysis of this, see the recent paper 
[iot). The Hodge star ★ is on the internal space and ([2|) 
gives rise to non zero supersymmetry breaking potential 
because of its lASD property. 

Once supersymmetry is broken, there would be non- 
trivial potential between the branes. The complete pic- 
ture can be worked out easily from the low energy type 
IIB supergravity lagrangian in the presence of D-brane 
sources in the following way (see also [ifj]): 
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non perturbative terms (3) 



where the dotted terms involve axio-dilaton plus other 
higher form interactions, Cq is a constant and = T ~ 
with T = — B. where is the gauge field on 
the D3 brane, Fi^ is the gauge field on the D7 brane, and 
{S, S') are the corresponding scalars on these branes. For 
more details see [l6i |. 

The first term in eqn. ([3]) when reduced to four di- 
mensions give rise to the kinetic terms for the complex 
and Kahler structure moduli. The second term is the 
potential term for the moduli: 
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GA*G 
Im T 



{\DW\^ ~ 3\W^\) (4) 



and therefore fixes (at least) all the complex structure 
moduli and the axio-dilaton (the 3|M^^| term goes away 
for the no-scale models). 

The second line in eqn. Q is related to the kinetic 
terms for the gauge fields on the D-branes and the in- 
fiaton S (if we ignore the scalar S"). The other scalar 
doublet X denote the hypermultiplet scalars with the cor- 
responding D-term and F-term potentials being given in 
the third and the fourth lines of eqn. ([3]). The potential 
AV, whose value will be specified later, denotes the ad- 
ditional attractive potential between the D3 and the D7 
brane. Finally, the non-perturbative potentials stabilize 
the remaining Kahler structure moduli. 
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The above picture is almost complete except for an 
important subtlety. The potential that fixes the complex 
structure moduli in eqn. ^ also gives a large mass to the 
inflaton field S because the Kahler potential K in eqn. 
(HI picks up a dependence on the inflaton ll|, |6l| . The 
potential for the inflaton then becomes very steep and 
consequently do not allow slow roll to happen in generic 
models. In the D3/D7 model this could be avoided (at 
least at the tree level) by allowing a shift symmetry, but 
this symmetry is broken by the AV term in ([3]). Thus 
a generic model of inflation from string theory is very 
complicated (see |63|, |6j| for some recent develop- 
ments) and since our aim here is to make the first step 
in understanding structure formation and reheating af- 
ter inflation, we will ignore the subtleties associated with 
moduli stabilization. 

Therefore, from the analysis presented above in eqn. 
([3]) we will only concentrate on the F-term and D-term 
potentials for the hypermultiplet scalars (which we will 
henceforth label as (/)+ and (f>- instead of the complex 
doublet x) and the inflaton S, and assume that the con- 
tributions from the background fluxes given in eqn. Q 
and the non-perturbative contributions fix most of the 
moduli. Note however that, as we mentioned earlier, 
there are other contributions to the inflaton potential 
coming from the Kahler potential X in (j4|) . For the sim- 
ple analysis that we present here, we will ignore this con- 
tribution. More details on this will appear elsewhere. 

Looking from the D3 point of view, the complex scalar 
S that parametrises the P^, and a set of complex scalars 
0+ and (jj- arising from the hypermultiplet-string be- 
tween the branes give rise to the following potential: 



y (5, (/.+ ,(/._) = 252|5|2(|^^|2 + |^_|2) 
+ 25' 
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C3)' (5) 



plus the additional AV term that we will add later. Here 
we have taken a more generic case by switching on and 
(3 as the FI terms to break susy. 

There is also another subtlety related to the issue of 
moduli stabilization that we should mention now. Due to 
background NS and RR fields, one might expect that all 
the complex structure moduli should be naturally fixed 
via a GVW potential [H] (see also [H, E^]). But our 
background is time-dependent, so the issue of fixing mod- 
uli should be addressed at any given time. Similarly, one 
might also expect that the supergravity no-scale struc- 
ture will be broken b y a llowing a gaugino condensate 
along the lines of [b^IsHI- These problems are not yet 
solved in a time dependent background, and an attempt 
to address these issues along the lines of [17| is under way 
(see also [H^ for another recent study of moduli correc- 
tions to D-term inflation). 

Finally, the issue of cosmic strings can be addressed 
by going to the full non-perturbative F-theory picture of 



the D3/D7 system (see 3 68| and the second paper of 
From F-theory it is easy to generate global symme- 
tries, and an SU{2) global symmetry can convert the cos- 
mic strings to semi-local strings of Achucarro-Vachaspati 
69I 70t type (see also 71 1). In terms of the D3/D7 sys- 
tem, this is simply a doubling of the number of D7 branes. 
These semi-local strings reduce the severity of the prob- 
lems generated by cosmic strings l7(| (they do not com- 
pletely eliminate the problems, as a recent study (t^ has 
shown) . 



THE VACUUM MANIFOLD OF THE MODEL 

From the previous section, we see that the potential 
which describes the dynamics of the complex inflaton 
field S (the separation of the two branes in the 4-5 plane) 
and of the two complex scalar fields </>+ and 0_ associated 
with the lowest energy modes of the strings stretching be- 
tween the branes is given by 

The vacuum manifold Ai is the set of field configura- 
tions which minimize the potential. The potential has 
been normalized such that the minimal value of the po- 
tential isV~Q. To be at the minimum of the potential, 
the inflaton field must have relaxed to 5 = 0. A special 
case which we will focus most of our attention on is the 
case (■+ = with (^3 7^ 0. Introducing radial and angular 
coordinates with each 6 field: 



p+e 



(6) 
(7) 



we see that the conditions to have V — become 



P+P- = 



The vacuum manifold M. is given by 

M : p+ = VCs , P- = . 



(8) 



(9) 



It is the manifold discussed in 17[. There are no 
non-compact flat directions. 

Another interesting special case occurs if = and 
C+ 7^ 0. In this case, setting the potential to zero yields 
the conditions 



(10) 



In this case, the vacuum manifold is the circle given by 

M : p+=P- = v/c72, ^+ = ■ (11) 
There are no non-compact flat directions. 
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In the general case there are no flat directions, either, 
and the vacuum manifold is given by 

2 2 

(12) 



P+ ^ 
^+ = 



In the following, we will focus mostly on the special 
case C+ = 0. 



SPECTRUM OF COSMOLOGICAL 
FLUCTUATIONS 

In this section we will review the calculation of the 
amplitude of the primary adiabatic perturbations. This 
is the mode which is usually considered in the literature 
and which is used to fix the model parameters. In Section 
7 we will turn to a discussion of the entropy mode. 

In the D3/D7 brane inflation scenario, inflation takes 
place while the branes are widely separated, i.e. while l^l 
is large. From ^ it follows immediately that for large 



values of the inflaton, the scalar fields 



and (i_ will 



vanish at the classical level (their quantum fluctuations 
will play an important role in the following section). 

Before turning on any supersymmetry breaking, and 
even after introducing the supersymmetry-breaking pa- 
rameters C+ and at the tree level, the configuration 
of branes is BPS, there is no force between the branes, 
and S will remain constant. One loop supersymmetry 
breaking terms will result in an additional one-loop con- 
tribution Ay to the potential 



Ay 
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(13) 



In the language of supersymmetric effective field theory, 
this term corresponds to D-term supersymmetry break- 
ing [t^. We are omitting the quartic contribution to 
the symmetry breaking potential which corresponds to F- 
term supersymmetry breaking. In the language of D3 /D7 
system, the above potential appears naturally from the 
supergravity solution of a D3-brane separated from a D7- 
brane by a radial distance r as : 



Ay cx arctan^ JF_ In 
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(14) 



where T- is related to the D7 brane gauge fluxes. There- 
fore the identification of ^3 with some components of the 
gauge fluxes on the D7 brane is consistent with expecta- 
tion. We also see that Sc in (|13p is related to the string 
scale A in (fT4)) . 

For small values of the string coupling constant g, this 
extra term in the potential will lead to a slow rolling of 



the inflaton field S towards 5 = 0. Initially, the values 
of (f)^ and (/>_ will vanish. Note that the presence of (fT5|) 
lifts the flat direction, preferring the value 0_ = 0. For 
a sufficiently small value of S, a tachyonic instability of 
the field will set in. The critical value \S\ = St below 
which this instability appears is given by 



(15) 



the condition coming from setting the second derivative 
of the potential (O to zero at =0. From direct open 
string calculation this is given by 



St 



Tra'arctan J-^ 



(16) 



Again we see that the identification of ^3 with some com- 
ponents of the gauge fluxes on the D7 brane is consistent 
with expectation. 

Since there is no instability in the (/>_ direction, 0_ 
will remain zero, again modulo quantum uncertainties 
which may become important during the later stages of 
preheating. We will therefore, for the moment, set 0_ = 
and consider the dynamics of the reduced configuration 
space consisting of the inflaton and . 

The inflationary scenario realized here is a speciflc case 
of hybrid inflation 7J| , and more precisely of supersym- 



metric hybrid inflation [73| (see also [75|). Whereas it is 
the field S which is slowly rolling and producing the den- 
sity fluctuations, it is the potential energy of the field 0+ 
which generates inflation. Fluctuations in general hybrid 
inflation models have been studied before (see, in partic- 
ular, [zilTi]). 

Using the standard theory of cosmologic al p erturba- 
tions in inflationary cosmology (see e.g. [771, IZSl for com- 
prehensive review articles and (79| for a pedagogical in- 
troduction) we obtain as amplitude 6h of the spectrum 
of scalar metric fluctuations 
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where rupi is the reduced Planck mass. In terms of the 
D3 /D7 model, this could be derived from the fluctuation 
of the background metric when the D3 brane moves to- 
wards the D7 brane. Now, making use of ([5]) and ([13]) we 
obtain (see also [l6'|) 
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5\/6 .gm3, 



(18) 



This is to be evaluated at the value of the inflaton S (we 
will from now on omit the absolute value sign for this 
field) when modes which are being observed today on 
large scales exit the Hubble radius during the period of 
inflation. For high scale inflation this happens ~ 50 
Hubble expansion times before the tachyonic instability 
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develops, at a value S — S^. Assuming that the acceler- 
ation of S during inflation is negligible, then S'jv is given 
by 



(19) 



where S is evaluated at = S't (we solve the scalar field 
equation of motion working backwards in time assuming 
constant S to obtain (fT9)) ). 
Provided that 



< 1, 



(20) 



the second term on the right hand side of pQ]) is smaller 
than the first. In this case, ([TS]) yields 
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3/2 



(21) 



On the other hand, if the inequality in (j20p is reversed, 
then the result becomes 
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(22) 



SCALAR FIELD DYNAMICS AND CONDITIONS 
FOR TACHYONIC RESONANCE 

The cosmological evolution starts with a separation be- 
tween the branes which is large on string scale. At large 
values of 15*1, the coupling between S and the string fields 
0-1- and 0_ constrains the latter to be zero, modulo clas- 
sical and quantum fluctuations. 

Neglecting the classical and quantum fluctuations of 
(j>+ and 0_ , the cosmological evolution is as follows: the 
inflaton field S rolls down the "valley" of its potential 
until a tachyonic instability for sets in. This occurs 
at the value \S\ = St given by (fT5|) . 

At that point, the evolution will take place in the com- 
plex two-dimensional field space of S and (j>^. There is 
no instability in 0_ direction, and hence 0_ will remain 
zero. At the background level, the term proportional to 
|(/)+0_| induces a linear confining potential for </)_. 

Let us now focus on the potential in (p^ direction for 
a fixed value of S. This potential has a negative effec- 
tive mass (i.e. a tachyonic instability) for small values 
of 1 0-1- 1 only. For larger values, the effective square mass 
is positive, leading to oscillations about the minimum of 
the potential. 

The phases of S and 0-|- do not play an important role 
in most of our considerations, and hence we will set them 
to zero in the following unless indicated otherwise. This 
will simplify the notation. We will briefly return to the 
issue of these phases in the concluding section. 



A new aspect of our work is the inclusion of the pos- 
sibility that the 0+ field has an offset from the symmet- 
ric p oint — at the end of the period of inflation 
[102| . Such an offset in any given Hubble patch will in 
general be produced by the fluctuations of the field on 



super-Hubble scales [81|, \8^. The offset will suppress 
the ability of small-scale quantum vacuum fluctuations 
of 0-1- to produce topological or non-topological defects 
on small length scales at the beginning of the resonance 
process [f03]. Such defect formation and the subsequent 
non-linear interaction of the defects was seen in previous 
work [4^ to be the most important aspect of reheating in 
hybrid-type inflation models like the one we are consider- 
ing. In our analytical study of reheating, we will assume 
that the effect of the small-scale vacuum fluctuations is 
sub-dominant compared to the effect of the offset. Even 
if the small-scale fluctuations dominate the reheating dy- 
namics, they will unlikely influence the evolution of the 
large-scale metric fluctuations which is the focus of most 
of our interest. 

The first question to ask concerning the effects of the 
offset is whether it can prevent the onset of the "tachy- 
onic" (or "negative coupling") [sl, [sl] resonance. For 
the sake of concreteness, we will assume that the offset 
is generated by the long wavelength fluctuations of (p+ 
which were generated during the period of slow-roll in- 
flation and check to see if they are sufficiently small such 
that, when S reaches the onset of the instability, 0-|_ will 
in fact be in the region of tach yon ic evolution. If this 
is the case then, as discussed in [46|, the energy transfer 
from the coherently evolving fields to matter fluctuations 
(the "preheating" phase of the reheating process) will be 
very rapid and complete before oscillations of 4>+ and S 
can set in. ff, on the other hand, the fluctuations of 4>+ 
are sufficiently large, then there will be no tachyonic pre- 
heating phase, and the energy transfer from the inflaton 
to the matter flelds will proceed by the parametric in- 
stability channel first discussed in [s^ (see also fsB]) and 
discussed comprehensively in [4fj . 

Before estimating the magnitude of the fluctuations of 
04-, we note that the slow rolling of the inflaton field S 
will prevent the onset of the instability of 0-^, until 



(23) 



Before this condition is satisfied, the evolution of the dy- 
namical system in the S — 0+ plane corresponds to a 
horizontal straight line (the solid curve in Figure f). 

We next find the region of phase space where the dy- 
namics is dominated by the instability in the (j)^ direc- 
tion. From the potentials ^ and (fT3|) (with 0_ — 0) it 
follows that the phase space curve where the two deriva- 
tives in (1231) are equal is given by the following relation: 



dV 


> 


dV 


90+ 




dS 



I 



2S 



C3 



g\3 

f67r2S' 



(24) 
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FIG. 1: Sketch of the evolution of the fields S and 0+ in 
configuration space. The solid line represents the evolution 
of the background fields. The cosmological evolution begins 
in the region of slow-roll (to the right and below the dashed 

curve). The slow-rolling ends when the dashed line is 
reached, after which a period of quasi-periodic oscillations 
about the ground state 5 = and 0+ = (^3 sets in. 



Above this curve, the dynamics is dominated by the 0+ 
instability. In the limit where (/>+ is small, we can simplify 
the potential such that the equation relating (f>^ and S 
has no quadratic or cubic powers of 0+ . If we also ignore 

1 /2 

the phase ip^ of and work in the limit S* ^ C3 (this 
approximation is good except at the beginning of the 
resonance) we get the following approximate relation: 



0+ 



9 



1 



(25) 



1/2 

The minimal value of cj)^ is taken on when S* ~ C3 , 
which yields a value of about p+ = 

7(16^2)^1/2^ The 

curve where the condition ([25| is satisfied is represented 
by the dashed curve in Figure 1. 

The condition for efficiency of tachyonic reheating is 



dpi 



< 0, 



(26) 



evaluated at the phase space point where the phase space 
trajectory of the background dynamics of our S — (j)+ sys- 
tem crosses the dashed curve given by (|25|). This would 
mean 



< 



Cs 2 



From here it is clear that the conditions 
correspond to 



P+ < 



3 ' 



(27) 
and (Ell) 

(28) 



which along the trajectory of Fig. 1 is automatically 
satisfied for small values of g if \S\ < St- 



We must now check whether the condition (pS]) is sat- 
isfied including the effects of fluctuations of 0+. Thus, 
we must calculate the magnitude of the root mean square 
fluctuations of 0+ . Since at the beginning of the cosmo- 
logical evolution, for large values of S, the valley of the 
potential is very steep, and since during inflation all pre- 
existing classical fluctuations are red-shifted, we expect 
the classical fluctuations to be negligible. Hence, we fo- 
cus on the quantum fluctuations. The quantum vacuum 
fluctuations averaged over a Hubble volume are given by 
integrating the Fourier space quantum vacuum fluctu- 
ations for all comoving modes which are super-Hubble 

m. 



< P+ >c 



(29) 



where the subscript Q refers to the fact that we are com- 
puting a quantum vacuum expectation value, and the 
He is the comoving Hubble constant. While the field S 
has a value much larger than St, the effective mass for 
the waterfall field 0+ will be large, and the vacuum fiuc- 
tuations small. However, if the effective mass remains 
smaller than the Hubble rate H for several Hubble ex- 
pansion times before the onset of the tachyonic instabil- 
ity, then the field fluctuations will approach the values 
obtained for a massless fleld. It can easily be verifled 
that the above condition for "effective masslessness" will 
be satisfied if (PO)) holds. In this case, one obtains 



< 



fHaW 

pI>Q- J d'k\Mk)\l 



,m— ' 



(30) 



where the subscript now indicates that we are taking the 
quantum fluctuations for a massless scalar field. Since 
in an expanding universe it is the field a{t)(j) which has 
canonical normalization, evaluating the right-hand side 
of (EH) gives 



< P+ >c 



2a?k 



(31) 



The above estimate is based on the same ideas which lead 
to stochastic inflation 'si'l. 

The Hubble expansion rate at the end of the slow-roll 
period is 



H ^ g 



pi 



(32) 



Plugging the value p+ ~ H with the above value of H 
into yields the following condition on the parame- 
ters of our model for which the reheating dynamics will 
go through a phase of tachyonic instability: 



C; 



1/2 



< 1. 



m. 



•pi 



(33) 
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We will initially choose the parameters of the D3/D7 
brane inflation model such that the direct adiabatic fluc- 
tuations produced during inflation, whose magnitude was 
discussed in the previous section, provide the primordial 
fluctuations observed today. It is now straightforward 
to check, making use of (I33|) and pT|) or ([22|) . whether 
tachyonic resonance preheating will occur. 

Let us consider some characteristic parameter values. 
Scales which are being observed today leave the Hubble 
radius about 50 e-foldings before the end of inflation. 
Hence, we take N = 50. In the range of values of g for 
which the inequality pp)) is violated, it follows from 
that for a COBE-normalized [sBl value of 5h — 10~^, we 
require 



.1/2 
^3 

mpi 
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-2.5 



(34) 



and that hence the condition for tachyonic reheating is 
trivially satisfied. For a value of g such as g = 10~^, the 
inequality ((20|) is violated. However, for sufficiently small 
values of g, there is a transition to the region in which 
pTjl rather than (|22|) applies. The conclusions concerning 
the validity of (j33[) . however, remain the same. 

To conclude this section, we have shown that the ini- 
tial stages of reheating will proceed, for reasonable pa- 
rameters of the D3/D7 brane system, via the tachyonic 
instability channel. However, once p+ has increased be- 
yond the value ([28| . the tachyonic decay will terminate. 
At this point, although most of the energy is contained in 
the nonlinear local fluctuations, a substantial fraction of 
the energy may remain in the background fields. To ob- 
tain an upper bound it takes on the time for most of the 
energy in the homogeneous fields to dissipate, we neglect 
the effects of the nonlinear fluctuations and analyze the 
subsequent dynamics which then proceeds by the para- 
metric resonance instability during the period when the 
background flelds oscillate about their minima. This is 
the topic of the following section. 

At this point, we must add some comments about the 
phase of S. Neglecting the phase of S means that the 
two branes (point particles in the X4 — X5 plane) will 
be approaching each other in a straight line. This re- 
quires fine-tuned initial conditions. Also, if these initial 
conditions were realized, then the D3 brane would most 
likely immediately be absorbed by the D7 brane when 

= for the first time. Thus, no oscillations of S could 
take place. In the following, we will be assuming that 
the impact parameter is sufficiently small such that the 
background field dynamics for S can be approximated 
by an oscillation of a real-valued field S. Oscillations of 
S about S — are interpreted as the D3 brane perform- 
ing approximately oscillatory motion about the D7 brane 
without zero impact parameter being reached. 



REHEATING IN THE D3/D7 BRANE 
INFLATION MODEL 



As discussed in the previous section, even though the 
initial stages of reheating in the D3/D7 brane inflation 
model take place by tachyonic resonance, this process 
shuts off early, leaving some of the initial energy in the 
form of oscillations of the background flelds. To obtain 
an upper bound on the time it takes to drain this residual 
energy from the background flelds, we proceed to study 
the phase when the background flelds S and cf)^ are oscil- 
lating about their respective minima. Under these condi- 
tions, preheating via a parametric resonance instability 
is expected to be the most efficient mechanism for energy 
transfer from the inflaton to matter field fluctuations. 

We recall that for <j>- — 0, the potential of our model 
reduces to that of supersymmetric hybrid inflation. Re- 
heating in hybrid inflation was flrst studied in j83 j. and 
then in greater depth in \M 



^ ^ ^ ^ where, in particular, 

the tachyonic instability channel was investigated lO'il. 
Some aspects of reheating in supersymmetric hybrid in- 
flation have been investigated in [88|. In the following, 
we take up this problem again and study it from a more 
analytical point of view. 

The strength of the preheating instability depends sig- 
nificantly on whether the resonance is "broad" or "nar- 
row", in the notation of [41]. Broad-band resonance is 
significantly more efficient. However, even narrow-band 
resonance can be much more efficient than the lowest 
order perturbative reheating channels first discussed in 
Below, we will find that, except possibly for a 
short period at the beginning of the phase of oscillations, 
we are not in the broad resonance region. Nevertheless, 
we will verify that the resonance is sufficiently efficient 
such that a substantial amount of the inflaton energy 
is transferred to matter fluctuations in one Hubble ex- 
pansion time. Following the philosophy first applied in 
[38|, we will neglect the expansion of space and use the 
efficiency condition just mentioned to argue that the ap- 
proximation of neglecting the spatial expansion is self- 
consistent. An improved analysis could be done by going 
to conformal time and rescaling the fields such that the 
Hubble damping terms in the equations of motion for the 
rescaled fields disappear. These equations have the form 
of Floquet-type equations and show exponential instabil- 
ities. In the absence of expansion of space, we obtain 
a particularly simple Floquet-type equation, namely the 
Mathieu equation. 

The analyses of preheating are based on studying 
the equations of motion for the fleld fluctuations semi- 
classically for a classical background dynamics. Before 
turning to the discussion of the fluctuation equations, we 
must specify the background evolution. As pointed out 
in Q , the background trajectories in the S — plane 
are typically not oscillatory along a line in the two dimen- 
sional configuration space. This introduces complications 



9 



compared to the usual studies of preheating. However, 
it was also shown in [sl] that the solutions are typically 
not too far from a preferred solution to the background 
equations where S and oscillate in phase (i.e. along 
a line in configuration space). Let us consider a partic- 
ular Hubble volume in which the fluctuation of has 
positive amplitude. In this case, we introduce a rescaled 
field 6 via 



/■1/2 
^3 



(35) 



Since the initial value of S at the end of the period of slow 
rolling is approximately St = y'Qi/2 whereas the ampli- 
tude of is almost vanishing, the approximate linear 
background trajectory during the resonance is character- 
ized by 



= -V2S. 



(36) 



We observe that both fields oscillate with the same fre- 
quency Lo which is given by 



25C: 



1/2 



(37) 



The potential of our problem is given by the sum of 
dS]) and (fTS]) . For small values of the string coupling 
constant g, the supersymmetry breaking term AV of (fT^ 
is negligible for small values of the inflaton compared to 
the first term ([5]). Hence, for 0_ = the relevant part of 
the potential becomes 



(38) 



where we recall that we are neglecting the phase of 0+ 
and thus treating (j)^ as a real variable. 

Now we are ready to derive the equations of motion for 
the field fiuctuations SS and S(p (to simplify the notation 
we are dropping the subscript -I- on S(p). We take the 
background field 



5*0 = -^Acos{Ljt) 
v2 

(l>o = Ca^^ - Acoa{ujt) , 



(39) 



where A is the amplitude of the background field oscilla- 
tions, and expand the potential to quadratic order in the 
field fluctuations defined by 



SS 



S — 5*0 



(40) 



Since we are expanding about a solution of the equa- 
tions of motion, the terms linear in the fluctuations in 
the action cancel by the background equations of motion. 
The quadratic terms in the potential are 



F(2)(55,, 



2g^(Pl6S^ + 8g^(f>QSo6S6(P (41) 



+2g'S'o 



Since the resulting equations of motion are linear, it is 
convenient to work in momentum space. The equations 
of motion for the fluctuation modes 6Sk and Sipk with 
comoving wavenumber k hence become 



dSk + {k'' + Ag^cl,l)dSk 
+8g'^(t)oSoS4'k 



(42) 



+8<720oS'oJ5fc = . (43) 

The amplitude A of the background held oscillations 
starts out of the order C,y^ but decreases as a conse- 
quence of the back-reaction of the fluctuations. As the 
fluctuations drain energy from the background, A de- 
creases. Since we are interested in whether preheating is 
efficient at draining most of the energy from the back- 
ground, we will in the following work in the approxi- 
mation that A? <C Cs ^-iid neglect terms quadratic in 
A. Inserting the background field expressions the 
equation for the inflaton fluctuation takes the form 



5Sk + {k' +^g''C3-8g''ACj^cos{Lot))5Sk 



\/2 



g'^ i^y^ Acos{ujt)5(l)k 



(44) 



Making use of the same approximations, the equation of 
motion for the matter field fluctuations 56 becomes 



5(j,k + {k" +Ag''C:i-l2g''AQ' cos{Lot))5(j,k 

= — ^g^Cl'^Acos{Lot)5Sk . 
v2 



(45) 



Fluctuations along the trajectory given by ([55)1 obey 

50fc = -V25Sk . (46) 



For these trajectories, both equations (|44|) and p5)l take 
the form 

5Sk + (fc' + Ag\^ - l&g-'ACJ^cos{ujt))6Sk = , (47) 

and can be put into the standard Mathieu equation form 

Xk + {Ak-2qcoB{2z))xk = Q (48) 

by introducing a rescaled time variable z via 

2z = cji, (49) 

and denoting the derivative with respect to z by a prime. 

By comparing and we can read off the values 
of q and Ak for the equation of motion for 5Sk'- 



* 1/2 ' 

S3 



(50) 



Ak = 4- 
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The condition for broad parametric resonance is 
Ak < 2q. 



(51) 



Hence, it follows from ([50)1 that, except right at the be- 
ginning of the reheating period whe n A"^ ~ Csj the broad 
resonance condition is not satisfied ^105| . 

We thus conclude that the resonance is of narrow-band 
type. This means 41, 91, Q^l that it will take place for 
values of k in resonance bands centered at half integer 
multiples of to. The lowest instability band is the widest. 
Its range is 



k e -[l-g/2,1 



9/2]. 



(52) 



We will now show that, in spite of the fact that the res- 
onance is narrow rather than broad, it is strong enough 
to drain a substantial fraction of the energy of the infla- 
ton within one Hubble expansion time, hence justifying 
neglecting the expansion of space in our analysis of the 
reheating equations. 

The condition for efficiency of the parametric reso- 
nance is 14 ll 



(53) 



where H is the Hubble expansion rate, evaluated at the 
beginning of the period of reheating. This condition 
comes from demanding that the exponential growth of 
the fluctuations induced by the parametric resonance in- 
stability is rapid compared to the Hubble expansion rate. 
The fluctuations grow exponentially, with the exponent 
being fi^t, where Hk is the so-called Floquet exponent 
given by 



fJ-k 



(54) 



The condition (|53p comes from demanding that fikSt > 
1 for the time interval St that a mode initially at the 
center of the resonance band will remain in the resonance 
band (recall that due to the Hubble expansion, modes are 
red-shifting with respect to the center of the resonance 
bands). 

Inserting the value of q from ([50|) and the value of H 
from ([32| , we find that the condition ((53|) for effectiveness 
of the parametric resonance instability becomes 



.1/2 ^ 72 



(55) 



which is easily satisfied for interesting values of Ca- 

We thus conclude that, in spite of the fact that pre- 
heating occurs in the narrow resonance region, it is suf- 
ficiently effective to drain a substantial fraction of the 
energy density of the background fields within a Hubble 
expansion time. 

The production of matter fiuctuations by the preheat- 
ing instability will drain energy from the background field 



oscillations and thus lead to a decrease in the amplitude 
A of these oscillations. The resonance will continue to 
be efficient until the condition ([55)) ceases to be satisfied. 
This will be the case once q decreases to the value 

qf = (-] , (56) 



which corresponds to a final amplitude of 



4 



1 C: 



1/2 



64 rupi 



<3 



(57) 



We immediately see that, provided ^ "^p/, 
stantial fraction / of the background energy density will 
transfer to fiuctuations during the time period when the 
resonance is efficient. This fraction is given by 



/ = 1- 



1 C: 



1/2 



(58) 



64 rripi 

We will complete this section with the discussion of 
two side issues. The first is the role of the mixing terms 
in the set of fluctuation equations (|42p . In principle, 
since the mass matrix for the fluctuations is symmetric, 
it could be diagonalized. Since the mixing terms are of 

1/2 

the order of A/(^ , the only effect they have would be 
to change the coefficients of the two q values slightly. 
Our conclusions about narrow versus broad resonance 
and about the efficiency of preheating are unchanged. 

The second issue concerns possible fluctuations in the 
0_ fleld produced during the preheating period. In the 
vacuum, 0_ = (see Section 3), and this is the value we 
will expand about. Note also that 0_ is real and posi- 
tive semi-definite. One can give a plausible explanation 
using duality arguments. Under a U-duality the string 
connecting the D3 and the D7 branes map to a wrapped 
D3 brane on a three cycle of a deformed conifold [93|. 
This is a massive black-hole in four dimensions and gives 
rise to a unique charged 4D hypermultiplet (941] . The 
expectation value of the scalars in this hypermultiplet 
gives the 4D black-hole configuration, and therefore it 
makes sense to take both 0± to be positive definite. This 
observation removes the apparent non-analyticity in the 
potential coming from the term proportional to 
As a consequence, every fluctuation about (/)_ = must 
be positive semi-definite in space and hence must have 
a non- vanishing spatial average. This average is subject 
to the linear confining potential mentioned above, and 
hence there is no possibility of a growth in fluctuations 
of 



ENTROPY AND INDUCED SECONDARY 
CURVATURE FLUCTUATIONS 

In this section we turn to the study of entropy fluc- 
tuations generated during the initial stages of reheating. 
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and to the calculation of the induced curvature fluctua- 
tions. The fact that entropy fluctuations induce a grow- 
ing curvature mode has been known from the early days 
of cosmological perturbation theory. An early applica- 
tion in the context of inflationary cosmology is to the 
computation of curvature fluctuations induced by axion 
perturbations in the early universe (see e.g. [95|). We 
will use the more recent formalism of Gordon et al. (96| 
for our analysis. 

We will focus on the effects of a single entropy mode 
and thus consider a two field system, the initial inflaton 
field S (which dominates the energy-momentum tensor) 
and the initial entropy field 0+ . The equation which de- 
scribes the growth of the induced curvature perturbation 
on super-Hubble scales is [96j 



n 



(59) 



where TZ is the curvature fluctuation in the comoving 
coordinate system, a is the effective adiabatic field, s 
is the effective entropy field, and Vg is the derivative of 
the potential with respect to the field s. The effective 
fields a and s are combinations of the inflaton field S 
and tachyonic field 0+ which depend on the background 
motion: 



a = S cos 6 + (hi sin 



and 



5s — S<j)+ cos 9 + SS sin 9 , 
where the angle 9 is given by 

cos - 



52 + 61 



(60) 



(61) 



(62) 



During the early stages of the tachyonic instability, the 
inflaton velocity is larger than the velocity of Hence, 
the angle 9 is approximately zero. 

The potential of our system is given by ([55)1 . Hence, 
the Hubble constant H at the onset of reheating is ap- 
proximately given by (|32p . 

During the tachyonic resonance period, the derivative 
Vb enteri ng in to our basic equation (|59p is approximately 
given by [106| 



-5^304 



(63) 



Since the Hubble friction is negligible, the equation of 
motion for the entropy field 0+ during this phase is 



which has exponentially growing solutions 
Mt) = </>+(0)e^-*, 



(64) 



(65) 



with Floquet exponent fip given by 



fJ-F 
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1/2 



(66) 



In the above, we set the time at the beginning of the 
instability to be t = 0. 

Neglecting the Hubble damping, the equation of mo- 
tion for the entropy fluctuation Ss ~ ^(f>-^- on super- 
Hubble scales becomes 



Ss 



Ss. 



(67) 



where the subscripts on the potential indicate which 
fields the derivative is taken with respect to. While the 
tachyonic instability condition (pS)) is satisfied, this equa- 
tion can be approximated by 



5s = g'^CsSs . 



(68) 



This shows that the entropy fluctuation increases with 
the same Floquet exponent as the background entropy 
field: 



Ss{t) = Ss{0)e' 



(69) 



The initial values of Ss and (f>^ are both given by the 
same large scale quantum fluctuations and will be taken 
to be of the same order of magnitude. 

Having determined the growth of the entropy mode, 
we can now integrate (I59p from the time t = when the 
instability starts to the time when the tachyonic instabil- 
ity condition breaks down. We denote the corresponding 
time hy t = tf. In the absence of back-reaction, the 
tachyonic resonance stops when (|28)) is satisfied. We de- 
note the result of the integration by ATZ. An approxi- 
mate evaluation of the integral (taking a to be constant) 
gives 



A7^(^/) 



1 



-(t>+{0)Ss{0)e 



2fj.pt f 



2mpi^F 

In the absence of back-reaction, then from 



.1/2 

^3 



(70) 



(71) 



and analogously for Ss{tf), we can solve for the duration 
t/ of the tachyonic instability. We see that the depen- 
dence of our result for A7^ on the initial value of the 
entropy mode drops out [l07l | and we obtain the result 



An 



AI2 




(72) 



Let us now turn to a brief discussion of back-reaction 
effects. The "waterfall" field 0+ has a dispersion a on 
microscopic scales which are due to its quantum vacuum 
fluctuations. By integrating up these quantum fluctua- 
tions of 0+ to its mass scale to, we obtain the following 
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result for the initial dispersion D{0) (that is the root 
mean square value of the field) at the time the tachyonic 
instability sets in 

D{0) ~ TO. (73) 

This dispersion then grows exponentially with the expo- 
nent set by the Floquet exponent fip, and after a time 

t. - ..In (^) (74) 

the dispersion will have grown to be comparable to the 
value of 0+ at the minimum of the potential. At this time 
(called the spinodal decomposition time), the field on 
small scales fragments into domains of typical size to~^. 
While in itself the formation of nonlinearities on small 
distance scales does not interfere with the linear growth 
of fiuctuations on cosmological scales, the nonlinearities 
can induce a positive contribution to the effective square 
mass of the tachyon field which shuts off the resonance. 

Whether back-reaction shuts off the resonance before 
the entropy fluctuation has had time to fully develop de- 
pends a lot on the values of the dispersion D(0) and of 
the quasi-homogeneous mode (/>+ on the scale k of the 
fluctuation. If the latter is calculated based on energet- 
ics, i.e. by setting the energy density in this mode to be 
comparable to the quantum vacuum energy density 

(j>lm^ - , (75) 

then the time scales tg and tf are comparable and our 
estimates hold. However, if 0+ is estimated by integrat- 
ing quantum vacuum fluctuations on length scales larger 
than k~^, and taking into account that the waterfall field 
is massive during most of the inflationary period and the 
amplitude of the fluctuations is hence redshiftcd, then 
ts is much smaller than t / , and hence the growth of the 
entropy fluctuations is cut off before they can become 
signiflcant. 

Let us for a moment assume that tf and ts are com- 
parable. In this case, our result ([7^ for the amplitude 
of the secondary curvature fluctuations induced by the 
entropy mode must be compared to the amplitude of the 
primary adiabatic mode. In the parameter regime where 
(f20|) is satisfied, i.e. 

Cg/^ > gmpi , (76) 
the primary adiabatic fluctuations are given by (see 

7^adiab ^ 5^'C3^'m;3. (77) 
Hence, in the parameter region given by 

grupi < Ca^^ < g^^^mpi (78) 

the secondary fluctuations dominate over the primary 
ones. 



When the condition ((20|) is not satisfied, the amplitude 
of the primary adiabatic modes is given by (see ([22| ) 

7^adiab ~ ViV^ , (79) 

and hence the secondary curvature fluctuations dominate 
over the primary ones unless 

VlVCg/' > mpi . (80) 

From the above results we can draw two main conclu- 
sions. Firstly, we see that the induced curvature fluctua- 
tion produced by the entropy mode remains in the linear 
regime. This is good news for the model. The second 
conclusion is that it is possible the secondary fluctua- 
tions dominate over the primary ones, thus necessitating 
a change in the model parameters in order to achieve 
agreement with the observed amplitude of fluctuations 
on large scales. 

DISCUSSION AND CONCLUSIONS 

We have studied reheating and structure formation 
in the D3/D7 brane inflation model of [l2|, with par- 
ticular emphasis on the tachyonic instability of super- 
Hubble scale entropy fluctuations. These entropy per- 
turbations induce a curvature fluctuations which we call 
"secondary" . We have found that under certain condi- 
tions these secondary fluctuations are larger than the pri- 
mordial adiabatic ones which have been considered in the 
past. They do, however, remain in the regime of appli- 
cability of linear perturbation theory. This would imply 
that the parameters of the model have to be changed 
compared to what is usually assumed in order to obtain 
agreement with the observed amplitude of the large-scale 
curvature fluctuations. 

Along the way, we have given an extensive analytical 
analysis of the reheating mechanism in the D3/D7 brane 
inflation model [lO^. In the low energy field theory limit, 
the dynamics of this system is a special case of supersym- 
metric hybrid infiation. We have taken into account the 
quantum fiuctuations in the "waterfall" field 0+. These 
play an important role for both the reheating dynam- 
ics and for the generation of entropy fiuctuations. For 
COBE-normalized values of the string theory parameters 
and for reasonable values of the string coupling constant 
we have found that, including the effects of the above 
mentioned quantum fluctuations, the initial stages of re- 
heating occur via a tachyonic instability. This instability, 
however, shuts off quite early, leaving some of the initial 
inflaton energy in the background fields. To obtain an 
upper bound on the time it takes to drain most of the 
residual energy from the background homogeneous fields, 
we have neglected the interactions of the nonlinear fluc- 
tuations produced during the initial phase of tachyonic 
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decay, and focused on the residual homogeneous field dy- 
namics. This later dynamics proceeds by the parametric 
resonance instability of [3^, [Hi [i^ In fact, except 
for at the onset of the reheating process, the system is in 
the narrow-band region of parameters. Nevertheless, the 
reheating is sufficiently efficient to convert a substantial 
fraction of the infiaton energy into matter fiuctuations. 

We note that efficient reheating in the D3/D7 brane 
inflation model is easier to achieve than in multi-throat 
inflation models since the matter flelds are directly cou- 
pled to the inflaton in the form of strings stretching be- 
tween the two branes whose separation constitutes the 
inflaton. 

In our work, we have neglected the issue of moduli sta- 
bilization. This is a very important caveat to our anal- 
ysis. Moduli stabilization in our model has been con- 



sidered in the first reference of 17| and more recently in 
97| . We plan to study the implications of the corrections 
to the potential induced by moduli stabilization on the 
dynamics of reheating in a foUowup paper. 

We have also simplified the dynamics of the back- 
ground fields. Namely, we have considered in-phase os- 
cillations of the two background fields S and , and we 
have neglected their phases. These phases provide ex- 
tra low mass entropy modes. It would be interesting to 
consider their excitation during reheating. 
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